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ABSTRACT 

Three observables of early- type galaxies - size (r e ), surface brightness (7 e ), and velocity dispersion 
(do) - form a tight planar correlation known as the fundamental plane (FP), which has provided great 
insights into the galaxy formation and the evolution processes. However, the FP has been found to be 
tilted against the simple virial expectation, prompting debates on its origin. In order to investigate the 
contribution of systematic stellar population variation to the FP tilt, we study here the FP relations 
of early-type galaxies in mid-infrared (MIR) which may represent the stellar mass well. We examined 
the wavelength dependence of the FP coefficients, a and b in log r e — a log (To + b log (J) e + c, using 
a sample of 56 early- type galaxies for which visible (V-band), near-infrared (if-band), and MIR 
(Spitzer IRAC, 3.6-8.0 /mi) data are available. We find that the coefficient a increases as a function 
of wavelength as da/dX — 0.11 ± 0.04 fim^ 1 , while the coefficient b reaches the closest to -1 at 3.6- 
5.8 /im. When applied to the visible FP coefficients derived from a larger sample of nearby early- type 
galaxies, we get the FP relation with (a, b) ~ (1.6-1.8,-0.9) at 3.6 /im. Our result suggests that the 
stellar population effect can explain more than half of the FP tilt, closing the gap between the virial 
expectation and the optical FP. The reduction in the FP tilt is reflected in the dynamical mass-to- 
light ratio, Mdyn/L, dependence on L which decreases toward 3.6-5.8 /im, suggesting that the MIR 
light better represents mass than the shorter wavelengths. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: formation — galaxies: fundamental 
parameters — galaxies: stellar content — galaxies: structure — infrared: galaxies 



1. INTRODUCTION 

In the search for correlations among physical parame- 
ters of early- type galaxies, it has been recognized that the 
effective radius (r e ), the effective mean surface bright- 
ness ((i) e )) and the central velocity dispersion (ao) 
form a planar relation (in logarithmic space) known 
as th e fundamental plane (here after FP; iDressler et al.l 
119871 : iDjorgovski fc David fl987t ). in the form of r e oc 
erg where a and b are foun d to be (a, b) ~ (1.2- 
1.5,-0.8) at visible w avelengths (|Jorgensen et alj 119961 : 
iBernardi et al.1 12003T ). Under the assumption of struc- 
tural homology and a constant mass-to-light ratio, the 
virial theorem implies that the FP coefficients should be 
(a, &)=(2,-l) - the so-called 'virial expectation'. The ob- 
served discrepancy, or tilt of the FP with respect to the 
virial expectation h as prompted many dis cussions to ex- 
plain its origin (see iD'Onofrio et~aT1 [20061 for a review of 
this field). 

One of the ideas is that the tilt is caused by the sys- 
tematic variation in the stellar population as a func- 
tion of physical pa rameters such as galaxy luminosity. 
iPahre et~a l. (1998b) investigated this effect by construct- 
ing the FP in X-band, which is supposedly a good tracer 
of the stellar mass less affecte d by age and dust. Mean- 
while, Scodeg gio et al.l (|1998f ) examined the wavelength 
dependence on the FP coefficients, and concluded that 
some of the tilt is caused by the stellar population man- 
ifested by the color-magnitude relation. These studies 
found that the stellar population effect exists, but it can 
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only partially explain the tilt of the FP. 

More recent investigations tackle the FP tilt prob- 
lem using new methods such as g ravitational lensing 
(iTreu et alj|2006t iBolton et al.ll2007h or dynamical mod - 
eling (jPadmanabhan et al l l2004HCappellari et al.l [2006) . 
Such studies suggest that the FP tilt nearly disappears 
when the FP is constructed in the mass-domain. Their 
conclusion is that the tilt must arise from the variatio n 
in the central mass-to-light ratio llRobertson et alf 2006) , 
but it is not clear wheth er the variation is d ominated by 
dark matter or by stars (|Bolton et al.ll2007h . 

In this Letter, we extend the FP analysis to wave- 
lengths beyond X-band to further investigate the effect 
of stellar population on the tilt. By doing so we aim 
to provide the FP that possibly better represents stellar 
mass (see § 5), and to improve the constraints on different 
scenarios for the FP tilt. 



2. THE SAMPLE 

Early-type galaxies were chosen from the sample of 
iPahrd (fl999h . which was used to study the FP of nearby 
early-type galaxies in visible and near-infrared (hereafter 
NIR) passbands. The sample has the velocity dispersion 
information necessary for constructing the FP. We then 
searched for mid-infrare d (hereafter M IR) 3 archival im- 
ages for galaxies in the iPahrd (|1999l ) sample. For the 
MIR data, we used the Spitzer Sp ace Telescope, Infr ared 
Array Camera (hereafter IRAC: [Fazio et ai1l2004h im- 
ages, covering four wavelength channels at 3.6, 4.5, 5.8, 
and 8.0 /im. The flux-calibrated, mosaiced IRAC images 
were retrieved from the Spitzer archive for these objects. 

The surface brightness fitting was performed for these 

3 We designate these wavelengths MIR to distinguish them from 
the K-b&nd. 
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TABLE 1 

Photometric parameters of the sample 



A 


^e,min 


fe,avg 


Te,max 




Mavg 


Mmax 


(Mm) 


(") 


(") 


(") 


(mag) 


(mag) 


(mag) 


0.55 


2.1 


20.1 


81.1 


-23.4 


-21.4 


-19.4 


2.2 


2.3 


14.8 


104.0 


-26.7 


-24.5 


-22.6 


3.(5 


2.0 


14.3 


65.3 


-25.8 


-23.6 


-21.8 


4.5 


2.2 


14.7 


80.5 


-25.1 


-23.0 


-21.2 


5.8 


1.2 


15.2 


90.0 


-25.1 


-22.7 


-21.0 


8.0 


1.1 


13.6 


86.8 


-24.2 


-22.3 


-20.5 



Note. — Effective radii and absolute magnitudes from Pahrc 1999 
(V- and X-band) and our Spitzcr IRAC analysis (3.6-8.0 /^m) are 
presented in minimum, average, and maximum values. 



matched galaxies, and the objects satisfying r e > 2" for 
three or more IRAC-bands were retained for the FP anal- 
ysis. We imposed this size limit in order to work with 
a sample with reliable r e values (see §3.1). After re- 
moving a few galaxies (NGC1275, NGC4824, NGC4478, 
NGC6166) that show peculiar light profiles (multiple 
source, close to a bright galaxy or stars), we finally iden- 
tified 56 galaxies with IRAC data in five clusters (A0426, 
A1656, A2199, A2634, and VIRGO) satisfying our selec- 
tion criteria. We present a brief summary of the photo- 
metric information in Table 1. The exposure times for 
the IRAC data range from 72 to 1000 sees. 

The above selection of the sam ple may introduc e a bia s 
in the derived FP coefficients (Scodegg io et al.l 1 1998h . 
However, such a bias would not affect our derivation of 
the wavelength dependence of the FP coefficients, since 
the multi- wavelength FP coefficients will be derived from 
the same galaxies for which the same bias would apply. 

3. ANALYSIS OF THE DATA 
3.1. Surface Brightness Fitting 

IRAF ELLIPSE was used to obtain surface brightness 
profiles of our IRAC sample galaxies. We restricted the 
fitting region to a > 2 pixels (along the semi-major axis) 
and discarded regions with S/N rms < 1. During the 
fit, we held the center, and fixed the ellipticities and 
the position angles of isophotes to those at the effec- 
tive radius in the 3.6 /im band. In addition, 3cr clip- 
ping was applied to reject outliers such as foreground 
stars. To subtract the background, we used the val- 
ues d etermined from the SExtractor (jBertin fe Arnoutsl 
1996). The adaptive background mesh sizes were varied 
between 16 to 96 pixels, and the best mesh was chosen 
to be the one which flattened the growth curve at the 
largest isophote (a ~3-6et e ). 

After the ELLIPSE photometry, we used the de Vau- 
couleurs r 1 ' 4 law to fit the observed surface brightness 
profiles measured along the semi- major axis 4 . The fit- 
ting procedure yields the effective radius (in ") r e = 
W (b/a) e a e where a e is the effective semi major axis and 
(b/a) e is the axis ratio of the isophote at this position. 
We tested the reliability of our fitting procedure using 
the simulated, PSF-convolved galaxies, and found that 
the surface brightness fitting gives unbiased, reliable re- 
sults when r e > 2". At the same time, we get the 
mean surface brightness within r e (in AB magnitudes) 
( M ) e - m 1/2 + 2.5 log (irr 2 e ) - 10 log (1 + z) - A x - K{z) 

4 We also tried the Sersic r 1 /™ law but found no difference in the 
FP coefficients; we therefore kept the r 1 / 4 methodology. 



where m 1 / 2 is the magnitude of the total flux within 
the effective isophote defined by a e and 6 e , while cosmo- 
logical dimming, galactic ex tinction (A\, using the for- 
mula of lLaureiis et al.lll996t and the extinction curve of 
iFitzpatrick fc Massal I2007T L and K-correction are taken 
into account. The K-correction is computed using the 
spectral energy distribution of a 13 G yr age, solar metal- 
licity, and 0.1 Gyr burst model from iBruzual fc Charlotl 
(2003), assuming the Salpter initial mass function. The 
last observable, a a is a kinematic parameter and is not 
expected to vary as a function of wavelength; we conse- 
quentl y use the sam e data used for the visible and NIR 
bands (jPahrelll999fl . 

In our analysis, angular sizes were converted into phys- 
ical length units for the FP construction by setting 
the distance to A1656 as 98.1 Mpc and calibrating the 
distances to individ ual clusters, utilizing the NIR FP 
(jPahre et a l. 1998a) as a distance ladder. 

3.2. Fitting of FP Coefficients 

We fitted the FP coefficients of the multi-waveband 
sample in the following manner using a variety of meth- 
ods: 

log r e = a log cr + b log(7) e + c, (1) 

where (/i) e and (I) e are related as (^} e oc —2.5 log(/) e . 
For the input r e and (I} e , we use our SB-fi t results for 
MIR (§ 3.1), and those listed in lPahrd (|1999h for V- and 
K-bands. We tried five different fitting methods: stan- 
dard least-squares fit, the inverse least-squares fit, the 
bisector of the two, the least-squares fit to the orthog- 
onal plane, and the least absolute deviation fit to the 
orthogonal plane. These methods are outlined below. 

It is natural to think of doing the s tandard least- 
squares fit of log r e (he reafter LSQ; iGuzman et al.1 
Il993t iBernardi et all l2003l ). but e arly FP work mainly 
took log a at the or dinate (|Dressler et al.l 119871 : 
iDjorgovski fc Davisl Il987t hereafter inverse LSQ) for 
their purposes. More recent work prefers the least- 
squares fitting of log r e by minimizing the variance or- 
thogonal to the FP pla ne (hereafter or t hogon al least- 
squares fit, or OLSQ; IBernardi et~aT1 I2003D or the 
least absolute deviations orthogonal to the plane (here- 
after orthogonal least absolute deviation fit, or OLAD; 
Pflrgensen etaT1ll996t IPahre et al.lll998af ). The orthogo- 
nal fitting has an advantage over other methods, reduc- 
ing the sys tematic error by treating the variables sym- 
metrically (jlsobe et al.l Il990f ) . However, the orthogonal 
methods yield larger measurement er rors than the LSQ 
methods, especially for small samples (|Isobe et al.lll990l) . 

Therefore, we also estimated the FP coefficients by tak- 
ing the bisector, or the plane equidistant from the planes 
obtained through the standard LSQ and inverse LSQ 
(hereafter the LSQ bisector). 1,000 Mo nte Carlo sam- 
plings of subsets of early- type galaxies in IBernardi et al] 
(2003) were performed to derive the FP coefficient errors 
on a sample of 50 early types to justify our approach. 
Through the sampling, we found the errors of the FP 
coefficients to be (cr a ,ab) — (0.14,0.06), best reproduced 
with the LSQ bisector method, while the other orthog- 
onal methods gave overestimated errors (>50% for the 
coefficient a). Aside from the error estimates, all three 
symmetrized method s reproduce the FP coefficient a of 
IBernardi et al.l (|2003l ) and the if-band early-type galaxy 
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TABLE 2 

Constructed Fundamental Planes at visible through MIR 



X(fim) 


a 


6 


c 


r 


(1) 


(2) 


(3) 


(4) 


(5) 


0.55 


1.23 ± 0.10 


-0.86 ± 0.04 


-9.16 ± 0.40 


0.96 


2.2 


1.42 ± 0.11 


-0.81 ± 0.05 


-8.20 ± 0.41 


0.95 


3.6 


1.55 ± 0.11 


-0.89 ± 0.04 


-9.89 ± 0.39 


0.96 


4.5 


1.47 ± 0.11 


-0.92 ± 0.04 


-10.16 ± 0.41 


0.96 


5.8 


1.57 ± 0.13 


-0.92 ± 0.05 


-10.55 ± 0.50 


0.95 


8.0 


1.55 ± 0.14 


-0.75 ± 0.05 


-9.30 ± 0.60 


0.93 


Note. 


— Fundamental planes for the 


sample of 56 galaxies with 


the r e > 


l" cut using the LSQ bisector method. Each column 


rcprc- 


sents (1) 


wavelength in /;m. 


(2)-(4) plane 


coefficients a, b, and 


c with 



associated uncertainties, and (5) the linear correlation coefficient. 



o- ^ 



6.2 6.4 6.6 6.6 7.0 7.2 
log u (kms~ 1 ) + 0.230</i> e 



7.0 7.2 7.4 7.6 7 
log t7 D {kms _, ) + 0.250</i> e 




8.0/xm 




6.6 6.8 7.0 7.2 7.4 7.6 
log tj a (kms" 1 ) + 0.236< J Li> e 



5.6 5.8 6.0 6.2 6.4 6.6 
log u a (kms -1 ) + 0.194</2> e 



Fig. 1. — Constructed fundamental planes projected in the di- 
rection of smallest scatter at 3.6, 4.5, 5.8, and 8.0 /im, respectively. 



sample of lPahre et al. (1998a) within 5% agreement. On 
the other hand, the standard and inverse LSQ methods 
are found to have about minus and plus 20 % system- 
atic biases in the coefficient a estimates respectively in 
comparison to the symmetrized methods. Given these 
results, we adopted the FP coefficients with the LSQ bi- 
sector method as our base result. 

4. RESULTS 

In Table 2, we list the FP coefficients with errors 
from 1,000 bootstrap resampling (unless obtained di- 
rectly from known error estimates, e.g., LSQ methods) 
derived from the LSQ bisector method, for wavelengths 
of 0.55-8.0 /im. We further plot the result of the FP 
fit in Figure 1, overlayed on the data points. We cau- 
tion readers to focus less on the absolute values of (a, b), 
but to focus more on the trend of the values with wave- 
lengths or methods (see discussions at the end of this 
section and § 2). Gathering the outcomes, we are led to 
the wavelength-dependent nature of the FP coefficients, 
with (a, b) values getting close to the virial expectation of 
(2,-1) as the wavele ngth increases. Such a tendency has 
been noted before (| Pahre et al.1 Il998bl : IScodeggio et al.1 
fl99cl . but our result indicates that it extends to 3.6 /im, 
and possibly beyond. When each cluster was analyzed 
separately, we also find the trend. 

The wavelength dependence of the FP coefficients is 
further presented in Figure 2, where they are plotted us- 
ing five different fitting methods (§3.2). The tilt of the 
FP is maximally reduced toward the virial expectation at 
IRAC-bands, and the thickness of the FP is maintained 
























"~~~^V virial 






— =4 plane - 



wavelength (/im) 

Fig. 2. — Wavelength dependence of coefficients a and b; dia- 
mond, triangle, square, cross, and plus sign symbols correspond to 
methodologies of inverse LSQ, OLSQ, LSQ bisector, OLAD, and 
ordinary LSQ respectively. Connected lines with spacing for identi- 
fication ar e from our catalog, while solitary symbols a r e from other 
literat u re IIDressler et aIJI1987fc IJgreensen et aLl fl996; Pah re et all 
I1998al : IBernardi et al.l 120051) of nearby samples with N > 100. 
Virial plane values assuming constant M/L are (a, fe)=(2, -1). 

thin for all but beyond 5.8 /im. In terms of the method- 
ology, we confirm the analysis of § 3.2 - the three sym- 
metrized fittings give coefficient a values that are fairly 
consistent with each other (considering the sample size). 

To quantify the wavelength dependence, we model the 
change of coefficient a as a linear function of wavelength 
by simultaneously fitting the OLSQ, OLAD, and the LSQ 
bisector results as follows: 



da/dX = 0.11 ±0.04 /^m" 



(2) 



from the visible to 3.6 /im (coefficient a behaves flat af- 
terward). This relation nicely explains the difference 
in coefficient o of .05 in the SDSS g*- to z*-bands 
(|Bernardi et al.l [2003f ). Meanwhile, for b, the tendency 
is not as linear as that for a, but has a local maxi- 
mum near the .ftf-band, approaches closest to -1 at the 
IRAC 3.6-5.8 /*m bands, and increases again at 8.0 /im. 
We attribute this behavior at 8.0 /im to the lower S/N, 
as well as the complexity in the 8.0 /im emission which 
can be dominated by the dust em ission in some cases 
(jBressan et alJl2006t IKo fc Imll2007t ). Indeed, the 8 .0 /im 
FP has the largest scatter among IRAC bands. The 
above result, together with the tendency of coefficient b 
from Table 2, implies that the increase in coefficient a and 
b are (Aa, Ab) ~ (0.34,-0.03) from V-band to 3.6 /im, 
and (Aa, Ab) ~ (0.15,-0.08) from if -band to 3.6 /im. If 
we use the FP coefficients from the references in Figure 
2 as the base values on which to apply equation (2), we 
obtain (a, 6) ~ (1.6-1.8,-0.9) at 3.6 /im, which is quite 
close to the virial expectation. The implication of this 
result is discussed in the next section. 

Note that our coefficient a in i f-band, derived from a 
subsample of 56 early- types from iPahre et alj (|1998af) is 
smaller than the value derived from their full sample of 
251 early-types by Aa = —0.11. The discrepancy should 
be mostly due to the limited sample size. More than 
half of our MIR galaxies belong to the Com a cluster (29 
objects ) , and the Coma cluster galaxies in I Pahre et al.l 
(1998a) show coefficient a in the if-band (a = 1.33) 
smaller than the total sa mple result by Aa = —0 .20, con- 
sistent with the results of iMobasher et all (|1999f ). Apart 
from the wavelength dependence, our results seem to be 
skewed to the FP of the Coma cluster. 

5. IMPLICATIONS ON THE ORIGIN OF THE FP TILT 

Recent studies suggest that the FP tilt originates 
mostly from a systematic variation in the mass-to-light 
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Fig. 3. — Left: The size vs luminosity test on the FP coefficient a. 
The solid line (labeled 'lum') is for the set of coefficients computed 
by exchanging the (I) e data of each wavelength with that from the 
K-band, while the dashed line (labeled 'size') is the similar result 
by exchanging the r e data. Right: The wavelength dependence of 
Mdyn/L on L as represented by the parameter (3 of M^ yn IL oc L' 3 
(see §5). 

ratio (jCappellari et alj 120061: iBolton et ai1l2007fh How- 
ever, the cause for the mass-to-light ratio variation is un- 
certain: it could be due to th e stellar population , or the 
central dark matter fraction (|Bolton et al.|[2007f) . Also, 
some studies suggest that the tilt is mostly explained by 
the non-homology related to the v ariation in the Sersi c 
index n among early- type galaxies (|Trujillo et al.l f2004). 
Here, we discuss the implication of our result on these 
issues. 

First, we investigated which one of the parameters - 
size or luminosity - dominates the observed change in 
the FP coefficients with increasing wavelength. This was 
done by deriving the FP coefficients from the MIR sam- 
ple by replacing (i) r e 's or (ii) (/} e 's, with those from the 
shorter wavelength data (in our case the K-band). The 
result is presented in Figure 3 (left), showing that the 
luminosity effect is the dominant factor up to 5.8 /im. In- 
terpretation at 8.0 /im is difficult due to low S/N and dust 
emission. Our result suggests that the stellar population 
effect is significant going from K-band to IRAC-bands. 

Next, we examined to what extent the stellar popu- 
lation plays a role in the FP tilt through the dynami- 
cal m ass-to- light ratio Md vn /L oc r e er§/L cx <jQ/(r e (I) e ) 
(e.g., iBernardi et al.1 12003) variation calculated from 
the FP coefficients. If Mdyn/L oc L 13 , then r e oc 

2A1+2/3) {1} - {1+mi+m _ ^ ^ of touji]lo et al | 

(|2004ft suggests (3 ~ 0.27 based on the visible FP. Our 
result is that the FP coefficient reaches a ~ 1.6-1.8 at 
3.6 /im (§4). In such a case, this relation gives (3 ~ 
0.06-0.13, which enables us to explain more than half 
of the tilt in the visible FP. Moreover, th e reduced tilt 
in th e mass plane (a^/p — off = 0.27: IBolton et al.l 
120071) is consistent with our Aa = 0.30 ± 0.11 from the 



Bernardi, M., et al. 2003, AJ, 125, 1866 

Bertin, E., k. Arnouts, S. 1996, A&AS, 117, 393 

Bolton, A. S., Buries, S., Treu, T., Koopmans, L. V. E., & 

Moustakas, L. A. 2007, ApJ, 665, L105 
Bressan, A., et al. 2006, ApJ, 639, L55 
Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 
Cappellari, M., et al. 2006, MNRAS, 366, 1126 
Djorgovski, S., & Davis, M. 1987, ApJ, 313, 59 
D'Onofrio, M., Valentinuzzi, T., Secco, L., Caimmi, R., & Bindoni, 

D. 2006, New Astronomy Review, 50, 447 
Dressier, A., Lynden-Bell, D., Burstein, D., Davies, R. L., Fabcr, 

S. M., Terlevich, R. J., & Wegner, G. 1987, ApJ, 313, 42 
Fazio, G. G., et al. 2004, ApJS, 154, 10 
Fitzpatrick, E. L., & Massa, D. 2007, ApJ, 663, 320 



I-band to the 3.6 /im in equation (2), advocating that 
the Mdyn/L variation is reduced by the regular light 
distributions in the MIR. As for the origin of the FP 
tilt, these results add another piec e of evidence against 
the s i gnificance of non-homology jP admanabha n et al.l 
120041: ICappellari eFaXl 120061 : IBolton et al.ll2007h . which 
predicts no change in the tilt with wavelength. 

We also derived the (3 parameter by directly fitting 
the Mdyn/L. Figure 3 (right) demonstrates that the ob- 
served dependence of Mdyn/L upon L decreases and be- 
comes flatter at IRAC-bands, just like the changes in (3 
derived from the FP coefficients. Combined with the fact 
that the change in the FP tilt with wavelength is dom- 
inated by the luminosity change, our Md yn / 'L-Gt result 
suggests that the rest-frame MIR luminosities at 3.6- 
5.8 /im better represent the stellar mass than the shorter 
wavele ngths, somewh at in agreement w ith iTemi et al.l 
(|2008l) , but not so with lLacev et all (|2007f ) . Among many 
possibilities, a proper combination of the metallicity and 
the age variation can possibly reproduce the observed 
trend, and we plan to investigate as future work, the 
physical origin of the Md yn /L - L relation as a function 
of wavelength. 

6. SUMMARY 

We studied the MIR fundamental plane relation of 
56 early-type galaxies and derived the wavelength de- 
pendence of the FP coefficients. When the FP is ex- 
pressed as r e oc ctq (I) b e) we found that the exponent on 
(j q , a, increases as a function of wavelength as da/dX = 
0.11 ± 0.04 /im" 1 , while b reaches closest to -1 without 
systematic variation with wavelength. When the FP co- 
efficients of previous studies are adopted as the starting 
point to calculate the MIR FP coefficients, we found that 
(a, b) ~ (1.6-1.8,-0.9) which is closer to the virial expec- 
tation, and that the change is dominated by the luminos- 
ity change. Together with the reduced dependence of the 
Mdyn/L on L at MIR wavelengths, our outcomes suggest 
that the variation in the stellar population is responsible 
for a significant portion of the FP tilt, and that the rest- 
frame MIR better probes the stellar mass of low redshift 
early- type galaxies than the shorter wavelengths. 

This study was supported by a grant (R01-2007-000- 
20336-0) from the Basic Research Program of the Korea 
Science and Engineering Foundation, and by the Seoul 
Science Fellowship (H J) . We thank the referee for useful 
comments, and Youngmin Seo and Soonyoung Min for 
algorithmic and technical advice in data analysis. 
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